Dissipation evidence for the quantum damped harmonic 

oscillator via pseudo-bosons 

F. BagarelloQ 



Abstract 

It is known that a self-adjoint, time-independent hamiltonian can be denned for the quan- 
tum damped harmonic oscillator. We show here that the two vacua naturally associated 
to this operator, when expressed in terms of pseudo-bosonic lowering and raising opera- 
tors, appear to be non square-integrable. This fact is interpreted as the evidence of the 
dissipation effect of the classical oscillator at a purely quantum level. 
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I Introduction 



In a recent paper, [I], we have analyzed the quantum damped harmonic oscillator within the 
context of pseudo-bosonic operators, which were introduced recently by Trifonov in [2] and 
further analyzed in a series of papers by us, [3]. The strategy adopted in pQ to quantize the 
system was that proposed in p[| following an idea given in [5] , idea which goes back to Bateman, 
[6], and which proved to be the most natural for our purposes. However, it is well known that 
Bateman's idea is not the only possibility of quantizing a generic dissipative system, and the 
damped quantum harmonic oscillator in particular: many different approaches exist and we 
refer to [7] for a detailed list of references. We just recall here that in quantum mechanics 
a dissipative system is quite often seen as an open system interacting with some reservoir, 
so that the energy of the system needs not be preserved, in general, and the unitary group 
of automorphisms which represents the dynamics for closed systems is replaced by a certain 
dynamical semigroup. 

Going back to our results in [I], we have shown that, even if pseudo-bosonic operators 
naturally appear in the analysis of the model, the vacua of the two annihilation operators of 
the model cannot belong to the Hilbert space "H = C 2 (M. 2 ) in which the system is assumed to 
live. Hence the general structure proposed in [3] apparently does not work. However, this is 
not the end of the story, since the details of our strategy are not at all uniquely determined. 
Thus it is natural to wonder whether different choices exist which allow us to construct the 
same functional settings as in [3]. In this paper we will continue this analysis, answering to this 
question. 

The paper is organized as follows: in the next section we introduce and discuss two- 
dimensional pseudo-bosons analyzing some of their mathematical properties. In Section III 
we introduce and analyze the quantum damped harmonic oscillator, while Section IV contains 
our conclusions. 

II The pseudo-bosonic settings 

Let "H be a given Hilbert space with scalar product (., .) and related norm ||.||. We introduce 
two pairs of operators, a, and bj, j = 1, 2, acting on "H and satisfying the following commutation 
rules 

[a,j,bj] = 1, and [ai,a 2 ] = [ai,b 2 ] = [a 2 A] = 0. (2.1) 
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Of course, they collapse to the OCR's for the independent modes if bj — aJj, j — 1,2. It is 
well known that aj and bj are unbounded operators, so they cannot be defined on all of %. 
Following [3], and writing D°°(X) := r\ p >oD(X p ) (the common domain of all the powers of the 
operator X), we consider the following: 

Assumption 1.— there exists a non-zero </?o,o £ H such that Oj</?o,o = 0, j = 1,2, and 

Assumption 2. there exists a non-zero \l/o,o £ H such that &j\I/o,o — 0, j = 1,2, and 
^o,o e J D 00 (a t 1 )n J D 00 (a t 2 ). 

Under these assumptions we can introduce the following vectors in %: 

= -i= 6? ^ ^o,o and # n ,, = -^(4)^ (4)% , n,l>0. (2.2) 
V^'<! V 7i! t! 

Let us now define the unbounded operators TVj := bjCij and 91-,- := AH = at&t, j = 1,2. It 
is possible to check that (/? ni / belongs to the domain of Nj, D(Nj), and ty n j G Z)(0Tj) , for all 
n, I > and for j = 1,2. Moreover, 

A^l^n.i = n^ n ,j, A^rM = ^n,«, ^l^n.i = ™^n,J, ^2^n,i = l^n,l- (2-3) 

Under the above assumptions, if we chose the normalization of \l/o,o and c£>o,o i n such a way 
that (^0,0,^0,0) = 1, we get 

($n,i, <fim,k) = 8n, m 5i,k, Vn, m,l,k> 0. (2.4) 

This means that the sets J 7 * = {^n,*; ^ > 0} and J-^ = {tp n .h n,l > 0} are biorthogonal and, 
because of this, the vectors of each set are linearly independent. We further assume that 

Assumption 3.— J 7 ^ and are complete in H. 

Let us now introduce the operators and Sty via their action respectively on the bases Fty 
and T^. 

S v ^ n ,k = Wn,ki Sv<fn,k = ^n,k, (2.5) 

for all n, k > 0. These imply that \l/ ni fc = (Sty Sy^n.fc and </? ni fc = (S<p Sty)(p n ,k, for all n, k > 0. 
Hence 

= S^ Sty = 1 =>■ Sty — Sp 1 . (2-6) 

In other words, both Sty and are invertible and one is the inverse of the other. Furthermore, 
we can also check that they are both positive, well defined and symmetric, [3]. Moreover, at 
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least formally, it is possible to write these operators in the bra-ket notation as 

oo oo 

S ¥> = ^2 \ < P*hk>< ( Pn,k\> S * = ^2 \^n,k><^n,k\- ( 2 ' 7 ) 

n,k=0 n,k=0 

These expressions are only formal, at this stage, since the series may not converge in the uniform 
topology and the operators S v and 5* could be unbounded. This aspect was exhaustively 
discussed in [3], together with many other features of pseudo-bosons which are not relevant for 
us here. 

It is interesting to observe that these two-dimensional pseudo-bosons give rise to interest- 
ing intertwining relations among non self-adjoint operators, see [8] and references therein. In 
particular it is easy to check that 

S* Nj = 9ljS* and Nj = % (2.8) 

j = 1,2. This is related to the fact that the eigenvalues of, say, N\ and 9Ti coincide and that 
their eigenvectors are related by the operators S v and S^. 

Ill Quantum damped harmonic oscillator 

In [TJ we have considered the quantum damped harmonic oscillator (QDHO) in connection 
with pseudo-bosons. Since the procedure is highly non unique, the negative results we have 
obtained in [1] only suggest the fact that something peculiar may happen. More explicitly, 
we have shown that the vacua of the pseudo-bosonic annihilation operators do not belong to 
£ 2 (R 2 ). This fact, however, leaves open the possibility that other choices are more appropriate 
than the ones considered in [1] , and in particular that the choice of the Hilbert space where the 
model should be considered was not the most appropriate. 

In this section we briefly review what we have done in [1] and then we propose different 
alternative approaches. All these proposals, however, give rise to the same conclusion: it 
seems to be impossible to have an Hilbert space in which the two vacua of the pseudo-bosonic 
annihilation operators live! We will comment on this result at the end of this section. 

IIL1 Working in C 2 (R 2 ) 

The original equation of motion of a classical damped oscillator, [5], mx + jx + kx = 0, is 
complemented by a second virtual equation, my — + ky = 0, and the classical lagrangian for 
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the system looks like L = mxy + ^(xy — xy) — kxy. This corresponds to a classical Hamiltonian 
H = p x x + p y y - L = ± (p x + 7|) (p y - 7|) + kxy, where p x = §§ and p y = §| are the 
conjugate momenta. The introduction of pseudo-bosons is based on two successive changes of 
variables and on a canonical quantization. First of all we introduce the new variables X\ and 
x 2 via x = -j^{xi + x 2 ), y = 75(^1 - x 2 ). Then L = \m (xj — x\) + 1 {x 2 ±i - Xxx 2 ) - \[x\ - x 2 2 ) 
and H = (y>i — 7"y) 2 + (P2 + 7^) 2 + f (^1 _ ^i)- The second change of variable is the 
following: 



P+ 

P- 
x + 

X- 



2mil 



■'mil 



Pi 
Pi 



X 2 , 



2u) + 



X\ + i 



1 



2w_ Xl 



X 2 , 

2m!)u + P 2 ' 
2mHZjr^2, 



(3.1) 



where we have introduced H = y ^ ^fc — and the two complex quantities u± = H ± z^. 

2 

In the rest of the paper we will assume that k > so that f2 is real. Up to now, we are still at 
a classical level, so that u + = p + = p_, x + = x_, and consequently, see below, H + = H_ 
and H = H. Hence H is a real Hamiltonian. Indeed, with these definitions, the Hamiltonian 
looks like the hamiltonian of a two-dimensional harmonic oscillator 

H = ~ (p 2 + + ulxl) + \ {pi + =: H + + H_ 

at least formally. 

At this stage we quantize canonically the system. Following jl], we require that the following 
commutators are satisfied: 

[x+,p+] = [x-,pJ\ = il, (3.2) 

the other commutators being trivial. We also have to require that p* + = p_ and that x\_ = x_, 
which are the quantum version of the compatibility conditions above. The pseudo-bosons now 
appear: 



(3.3) 
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and indeed we have [a+,6+] = [a_,6_] = 1. Notice also that b + = a_ and 6_ = a + . Moreover 
H can be written in term of the operators N± = b±a± as H = u + N + + u.N- + ^±±^ 1. So 
the hamiltonian of the QDHO is quite simply written in terms of pseudo-bosonic operators. 
In [1] we have used the following representation of the operators in (13.2]) : 

x + = TTTsT (pPv + t*) ' 
< x - = ¥J^f (FPy + 6x), (3 4) 

p + = Tp x + 5y, 
^ p~ = Tp x + Sy, 

for all fixed choices of V and 8 such that T S ^ 5T. Here sc, y, p x and p y are pairwise conjugate 
self-adjoint operators: [x,p x ] = [y,p y ] = £1. Notice that these operators also satisfy the 
compatibility conditions p\ = p- and x + = x~. Hence it is natural to represent x and y as the 
standard multiplication operators and p x and p y as — % -J^ = —i d x and ~^-§^ — —id y . Then we 
get 




iady 
iad„ 



— d x + iad v 



— d x + iad v 



(3.5) 



where, to simplify the notation, we have introduced a = T , T r c ^ and 3 = T , T <5 -= . Notice that, 
r J r<5— 8r 1 r<5-<5r 

since these operators satisfy the pseudo-bosonic commutation rules, the coefficients in (13. 5p 
satisfy the equalities a 5 = j3 Y and 3 T — a 5 = 1 . 



Due to the fact that b A 



a_ and b. 



A, 



Assumptions 1 and 2 of Section II collapse 



and we just have to look for a single square-integrable function ip 0) o(x,y) such that, first of 
all, a + cpo >0 (x, y) = a-(p 0j0 (x,y) = 0. It is possible to check that a solution of a + (p 0)0 (x,y) = 
a^Lpo :0 (x,y) = is the following: 



<po,o{ x >y) = N f ex P 



2T 
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x 



2au + 



(3.6) 



where it is required that the equality — = 
both 0+^0,0(^5 y) = and a_ip fi(x, y) = 0. Of course, we would like (po,o( x , y) to be square 



= holds, since in this way ^0,0(^5 y) satisfies 



integrable. It is easy to check that 



I3LO+ 



and 



are real. This follows from the above relation 
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between u + and Hence <po,o{x,y) G C 2 



only if 



^ 4 



> and if, at the same time, 



a u;_|_ 



< 0. However, using the previous relations between the coefficients, it follows that 



6£ 



> 0. 



Hence (fo,o( x ,y) cannot belong to C 2 (R 2 ). We conclude that Assumptions 1 and 2 of Section 
II are not satisfied, so that apparently there is no possibility of constructing two biorthogonal 
bases of C 2 (R 2 ) out of the QDHO. 



III. 2 Changing Hilbert space 

This result does not exclude, however, that a solution of Assumptions 1 and 2 could be found 
in a different C 2 space, for instance in a space with a suitable weight. A simple-minded idea 
would be to replace C 2 (R 2 ) with, for instance, Hi := £ 2 (M 2 , e~ Clx2 ~ C2y2 dxdy), where c\ and c 2 
should be two positive constants chosen in such a way that the wave-function ip ^(x, y) in (]3.6p 
does belong to Hi- This approach, however, has an immediate drawback: in Hi the adjoint 
of the operator is different from the one in £ 2 (R 2 ), and this should be taken into account to 
produce a consistent model. This is what we will do in this section. In particular we will show 
that what we expected is reasonable but false! 

In Hi the scalar product is clearly defined as follows: (/, g) 1 = L dx L dy f{x) g{x) e ~ Clx2 ~ C2V 
The adjoint X* of the operator X in Hi, is defined by the equation (Xf,g) 1 = (f,X*g) v for 
all f,g& Hi belonging respectively to the domain of X and X* . We refer to [9] for further and 
rigorous reading on the definition of the adjoint for unbounded operators. It is easy to check 
that d* = —d x + 2cix and d* = —d y + 2c 2 y, which return the adjoint in £ 2 (R 2 ) if ci = c 2 = 0. 
Taking as our starting point the operators a± and b± defined in f !3.5p . we can compute their 
new adjoints (i.e. their adjoints in Hi), which can be written as 

a* + = 6_ + 
a*_ = b + + ^J2oo + 
b* + = a_ + y/2uT 
b*_ = a + + 

It is clear that these again reduce to the adjoints in £ 2 (R 2 ) when c\ = c 2 = 0. It is also clear 
that, but for this case, the compatibility conditions required above are not satisfied: a?j_ ^ b T . 
Nevertheless, if we carry on our analysis, we can still look for the solutions of the differential 



cix + ic 2 ya) 
cix ^ + ic 2 yaj 
-cix + ic 2 ya 
—c\x + ic 2 ya 



(3.7) 
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equations a + (po t o(x,y) = a_Lp 0)0 (x,y) = and b^^fo^(x,y) = b*_^o,o{x,y) = 0. The solution 
fo,o(x, y) is, clearly, exactly the one in ( 13.61) . with the same condition on the ratio ^ as before. 
The wave- function <p $(x,y) belongs to Hi if the following inequalities are satisfied: 

ci + ^± > 0, c 2 - — > 0. (3.8) 

Notice also that ipo t0 (x,y) is eigenvector of H with eigenvalue \ (oj + + ujJ). This may appear 
as an improvement with respect to the result of the previous section, since non trivial choices 
of Ci and c 2 for which (13.81) are satisfied do exist. For any such choice <pofl(x, y) belongs to Hi 
and it also belongs to the domain of all the powers of 6_ and b+, so that the wave- functions 
Vn + ,n_ (x, y) can be defined as in (I2.2p . Let us now look for the function v&o.o^j y)- Due to (13. 7p 
we get 

*o,o(z,y) = N * ex P (- ^rrr ^ + tt— — 2/ 2 }exp {cix 2 + c 2 y 2 } , (3.9) 
[ 1 1 la oj + J 

which belongs to "Hi if 

^-±-Ci>0, c 2 + <0. 3.10 

^o,o(%,y) is eigenvector of i/*, which is different from W = H, with eigenvalue | (cu + + o;_) 
It is not hard to check, now, that these conditions are not compatible with those in ( I3.8p : in 
other words it is not possible to fix Ci and c 2 in such a way both ( 13. 8 j) and ( 13 . 1 0[) are satisfied. 
This means that our original simple-minded idea that adding a weight in the scalar product of 
the Hilbert space should regularize the situation does not work as expected. More explicitly, if 
Ci and c 2 satisfy ( I3.8p . then ^ 0t0 (x,y) does not satisfy Assumption 2. Viceversa, if they satisfy 
(I3.10p . then ^0,0(^5 y) does not satisfy Assumption 1. In both cases, therefore, only a single set 
of functions in "Hi can be constructed, which is (most likely) a basis of Hi itself. 

III. 3 And now, what? 

Changing Hilbert space has proven an interesting exercise but it doesn't give positive results. 
For this reason, in this section, we stay in the original Hilbert space and we consider four 
operators a± and b± satisfying b± = a) T and the following commutation rules, for convenience 
given only in terms of a± and of their adjoints: 

[a + ,ol] = 1, [a+,a_] = [a + ,a+] = [o_,a_] = 0. (3-H) 

A rather general expression for a±, which extends the one in (13. 5p . is the following 

a + = a x x + a y y + (3 x d x + (3 y d y , 
a- = -y x x + 7 y y + rj x d x + %d y , 
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where a x , a y , . . ., r) y are complex constants to be fixed. Due to the form of the operators, it is 
natural to look for solutions of equations a±fo ! o(x, y) = in the form fo,o(x, y) = N e~ kix2 ~ k2y2 , 
for some positive k\ and k^. This choice would guarantee that fo,o{x,y) belongs not only to 
£ 2 (M 2 ), but also to D°°(b + ) n D°°(b_). However, it is easy to check that, with these choices, 
no solution exist: simple algebraic manipulations show that /3 X — j3 y = r] x = rj y = and, 
consequently, that also a x = a y = j x = ^ y = 0, which is not possible. Hence, if we want 
to keep the form of the operators as in (13.121) . we have to change the expression of fo,o(x,y). 
The simplest extension is fo t o(x,y) = N e~ klx2 ~ k2y2 ~ k3Xy , where k 3 is a third real or complex 
constant to be fixed and which does not have to destroy the square-integrability of /o,o- It is 
possible to check that, if a solution does exist, it must correspond to non real values of the /3's 
and of the 77's. Indeed it is again a matter of simple algebraic computation to check that, if 
Pxi Py, Vx and rjy are assumed to be real, then no solution of the system does exist. It is also 
possible to check that, if k 3 = 2\fk\ &2, no solution exist both for real and for complex values of 
the coefficients in (I3.12p . In other words: in all the situations which we have completely under 
control the solution simply does not exist, in the sense that we may have square integrability 
in, say, x but not in y or vice-versa. Or, yet, we can have square integrability of (pofl(x,y) 
in both x and y but not of the second wave-function ty 0t0 (x,y). It is also clear that changing 
the definitions in (13.121) is dangerous since it can destroy the compatibility conditions or the 
commutation rules, or both. Analogously, changing the expression for fo,o( x ,y) could easily 
destroy the validity of Assumption 1: our choices seems to be the only possible ones! 

At a first sight these results may appear unpleasant. On the other hand, in our opinion they 
are very much physically motivated. Indeed, let us recall our starting point: our real system, 
a damped harmonic oscillator, was coupled to a virtual, somehow forced, harmonic oscillator. 
Hence the energy lost, at a classical level, from the first oscillator was transferred to the second, 
in such a way that the full system stays conservative. This classical flux of energy produces, 
at a quantum level, the impossibility of having both </?o,o and \l/o,o in a single Hilbert space: 
for instance, if 990,0 can be normalized, \I/o,o can not: its norm increases too much! And this 
result (or the related ones) seems to be independent of several possible way out that one may 
imagine: if we change Hilbert space to allow more functions in, since we also have to change 
the adjoint, we go out of the space anyway. Or we can look for other, and possibly more 
convenient, representations of the pseudo-bosonic operators. But, once again, requiring that 
their vacua are square-integrable produces a set of algebraic equations with no solution at all. 
Needless to say, our choices do not cover all the possible situations, so that, in principle, other 
rather special solutions reflecting the construction in Section II could be found. However, we 



9 



claim that these solution cannot exist, exactly because of the non conservative nature of our 
oscillator. 

IV Conclusions 

In this paper we have considered in more details the QDHO already discussed in [I]. We have 
shown that working in a different Hilbert spaces or changing the representation of the pseudo- 
bosonic ladder operators do not change the main results of our previous analysis, i.e. that the 
formal eigenstates of the number operators are not all square-integrable. This impossibility has 
been interpreted here as intrinsically related to the non conservative nature of the quantum 
system under analysis. 
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